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Abstract. Random matrix theory can be used to describe the transport properties 
of a chaotic quantum dot coupled to leads. In such a description, two approaches 
have been taken in the literature, considering either the Hamiltonian of the dot or its 
scattering matrix as the fundamental random quantity of the theory. In this paper, we 
calculate the first four moments of the distribution of the scattering matrix of a chaotic 
quantum dot with a time-dependent potential, thus establishing the foundations of a 
"random scattering matrix approach" for time-dependent scattering. We consider the 
limit that the number of channels N coupling the quantum dot the reservoirs is large. 
In that limit, the scattering matrix distribution is almost Gaussian, with small non- 
Gaussian corrections. Our results reproduce and unify results for conductance and 
pumped current previously obtained in the Hamiltonian approach. We also discuss an 
application to current noise. 
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1. Introduction 

From a statistical point of view, energy levels and wavefunctions in semiconductor 
quantum dots and metal grains, or eigenfrequencies and eigenmodes of microwave 
cavities, share a remarkable universality. With proper normalization, correlation 
functions of energy levels or wavefunctions for an ensemble of macroscopically equivalent, 
but microscopically distinct samples depend on the fundamental symmetries of the 
sample only; they do not depend on sample shape or volume, or on the impurity 
concentration. The same universality appears for correlators of eigenvalues and 
eigenfunctions of large matrices with randomly chosen elements |^, |^. Originally, 
such "random matrices" were introduced by Wigner and Dyson to describe the universal 
features of spectral correlations in heavy nuclei 0, Theoretical predictions from 
random matrix theory have been verified in experiments on semiconductor quantum dots 
and chaotic microwave cavities, and with the help of numerical simulations p, |^, p, 
for the case of a disordered quantum dot, the validity of random matrix theory has been 



proven by field-theoretic methods [10 



Open samples, such as semiconductor quantum dots coupled to source and drain 
reservoirs by means of ballistic point contacts or microwave cavities coupled to ideal 
waveguides, do not have well-resolved energy levels or wavefunctions. They are 
characterized by means of a continuous density of states and by their transport 
properties, such as conductance or shot noise power. Within random-matrix theory. 



two approaches have been taken to describe open samples ||rT|. In both approaches, 
transport properties are described in terms of the sample's scattering matrix S. The 
first approach is the "random Hamiltonian approach" . In this approach, the scattering 
matrix is expressed in terms of a random Hermitean matrix H, which represents the 
Hamiltonian of the closed sample. Averages or fluctuations of transport properties are 
then calculated in terms of the known statistical distribution of the random matrix 
H |1T2| , In the second approach, the "random scattering matrix approach", the 



scattering matrix S itself is considered the fundamental random quantity. It is taken 
from Dyson's "circular ensemble" of uniformly distributed random unitary matrices 



|TJ, |T5|, |T6|, |T^, |T8[, or a generalization known as the "Poisson kernel" P, |T9|, Both 
approaches were shown to be equivalent |jT3|, 0. Hence, in the end, which method to 
use is a matter of taste. 

Recently, there has been interest in transport through chaotic quantum dots with a 
time-dependent Hamiltonian. Switkes et al. fabricated a "quantum electron pump" 
consisting of a chaotic quantum dot of which the shape could be changed by two 



independent parameters I^Tj. Periodic variation of the shape then causes current flow 
through the quantum dot, hence the name "electron pump". Motivated by theoretical 
predictions of Vavilov and Aleiner [^, Huibers et al. looked at the effect of 
microwave radiation on the quantum interference corrections to the conductance of 
a quantum dot The presence of a time-dependent potential will cause the ratio 



of universal conductance fluctuations with and without time-reversal symmetry to be 
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less than two if the typical frequency of the fluctuations is of the order of the electron 
escape rate from the quantum dot p3| , p5| , p6| , p7| . (By the Dyson-Mehta theorem |28 



the ratio is two in the absence of a time-dependent perturbation [11 



A scattering matrix formalism to describe time- dependent transport was developed 
by Biittiker and coworkers |2^, 0, ^ The scattering matrix formalism for 
time-dependent scattering is more complicated than the formalism for time-independent 
scattering, since energy is no longer conserved upon scattering from a cavity or quantum 
dot with a time-dependent potential. In the adiabatic limit, when the frequency u of the 
time-dependent variations is small compared to the escape rate from the quantum dot, 
the theory can be formulated in terms of the scattering matrix S and its derivative to 
energy . However, a theory that describes arbitrary frequencies u must be formulated 
in terms of a scattering matrix S{e,e') that depends on two energy arguments, or. 



equivalently, a matrix S{t,t') depending on two time arguments P2 



Random matrix theory can be used to describe the statistics of time- dependent 
transport if the time dependence is slow on the scale of the time Terg needed for ergodic 



exploration of the quantum dot. In several recent papers |22, 23, 25, 26, 27, 34, the 



calculation of time-dependent transport properties for an ensemble of chaotic quantum 
dots was done using a variation of the "random Hamiltonian approach": the time- 
dependent scattering matrix S(t,t') is first expressed in terms of a time-dependent 
Hermitean matrix H{t), which is the sum of a time-independent random matrix and 
a time-dependent matrix which does not need to be random; the ensemble average is 
then calculated by integrating H over the appropriate distribution of random matrices. 
It is the purpose of this paper to develop a "random scattering matrix approach" for 
time-dependent scattering, using the distribution of the scattering matrix S{t,t'), not 
the Hamiltonian H{t), as the starting point for further calculations. 

For time-independent scattering, the distribution of the elements of the scattering 
matrix is given by the circular ensembles from random matrix theory (or, for a quantum 
dot with non-ideal leads, by the Poisson kernel). In the limit that the dimension of 
the scattering matrix becomes large, the scattering matrix distribution can be well 
approximated by a Gaussian, whereas non-Gaussian correlations can be accounted for 
in a systematic expansion in 1/A^ ||3^. For the calculation of transport properties 
(conductance, shot noise power), the Gaussian approximation is usually sufficient; 
knowledge of the underlying "full" scattering matrix distribution is not required. Here, 
we take a similar approach for time-dependent transport. We show that, for large A^, 
elements of the scattering matrix S{t,t') are almost Gaussian random numbers, for 
which non-Gaussian correlations can be taken into account by means of a systematic 
expansion in 1/A^. We calculate the second moment of the distribution and the leading 
non-Gaussian correction. 

This paper is organized as follows: In section Q we review the scattering matrix 
approach for time-independent scattering. The case of time-dependent scattering is 
considered in section ^. Applications are discussed in section ^. Details of the calculation 
and an extension to the case of quantum dots with nonideal contacts can be found in 
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the appendices. The second moment of the scattering matrix distribution calculated 



here was used in reference |^ to compute the shot noise power of a quantum electron 
pump. 

2. Time-independent scattering 

We first summarize important facts about the distribution of the scattering matrix S 
for time-independent scattering. 

For large matrix size A^, the scattering matrix elements Sij have a Gaussian 
distribution with small non-Gaussian correlations. Mathematically, this is a consequence 
of the fact that S is distributed according to the circular ensemble from random matrix 



theory or, for a quantum dot with non-ideal leads, the Poisson kernel |3g]. In a 
semiclassical picture, the Gaussian distribution of the scattering matrix elements follows 
from the central limit theorem, when Sij is written as a sum over many paths, where 



the contribution of each path contains a random phase factor [|TJ, |T5|, [T^, |3^, |3^ . The 
small non-Gaussian corrections follow because the full scattering matrix satisfies the 
constraint of unitarity, which is not imposed in the semiclassical formulation.||] 

The Gaussian part of the distribution is characterized by the first two moments. In 
the main text, we focus on the case of a quantum dot coupled to the outside world via 
ideal leads. In this case, the first moment vanishes, 

{S.,) = 0. (1) 

The case of nonideal leads, for which {Sij) ^ 0, is discussed in appendix A. The second 
moment of the scattering matrix distribution depends on the presence or absence of 
time-reversal symmetry (TRS), 

U^iJ-M _ / c C* \ _ ^ X. J i^ik^jl + ^il^Jk) ^Ith TRS, , . 

W, - - ]^ X j s^j^^ without TRS, 

up to corrections of relative order in the presence of time-reversal symmetry. All 
averages involving unequal powers of S and S* vanish. Equation (|]) is for spinless 
particles or for electrons with spin in the absence of spin-orbit coupling. (In the latter 
case, the scattering matrix has dimension 2N and is of the form S ® I2, where I2 
is the 2x2 unit matrix in spin space and the N x N matrix S describes scattering 
between orbital scattering channels.) We do not consider the case of broken spin-rotation 
symmetry, when iS is a random matrix of quaternions 

Non-Gaussian correlations of scattering matrices are of relative order 1/A^ or less. 



The leading non-Gaussian correlations are described by the cumulant |^ 

JJ/hjl,i2j2',k\ll,k2l2 / C C C* C*\ /C C*\/C C*\ 

^^^2 — \<-'lljl<-'i2j2<-'kih<-'k2l2/ ~ \'-'hjl<-'kih/Wi2j2'~'k2l2/ 

I Averages or correlation functions of certain transport properties which, at first sight, would require 
knowledge of the fourth moment of the scattering matrix distribution, can be formulated in terms of the 
second moment only, using unitarity of the scattering matrix. This way, the average and fluctuations 
of the conductance of a chaotic quantum dot have been calculated using the semiclassical approach, 
see, e.g., references [40, 



Scattering matrix ensemble for time- dependent transport 



5 



- {<Shj^<Sl^i2) {<Si2j2<Skih) ■ (3) 
In the absence of time- reversal symmetry and for ^ 1, the cumulant W2 is given by 

In the presence of time-reversal symmetry, W2 is found by the addition of 14 more terms 
to equation corresponding to the permutations 12 ^ j2, ki ^ li, and k2 ^ h, 

W2 = i^hki^jil2^i2k2^j2h + ^hk2^hh^i2kx^j2l2 + ^ilkl^jll2^j2k2^i2h + ^hk2^jlh^j2kl^i2l2 

~l" ^iih^jil2^i2k2^j2ki ~l" ^iik2^jiki^i2li^j2l2 ~^ ^iih^ jih^ j2k2^i2ki 
~\~ ^iik2^ j\ki^ j2h^i2h ~^ ^iiki^jik2^i2h^j2h ~^ ^iil2^jih^i2ki^j2k2 
+ ^iiki^jik2^j2l2^i2h ~^ ^iil2^jih^j2ki^i2k2 ~^ ^iih^jlk2^i2h^j2ki 

+ ^iih^ jiki^i2li^ j2k2 + ^iili^jik2^j2l2^i2ki + ^iil2^jikiSj2h^i2k2) ■ (5) 

We refer to reference for higher-order cumulants and finite- corrections to Wi and 
W2. 

Although equations (H) and @ do not specify the full scattering matrix distribution 
— for that one would need to know all cumulants — , they are sufficient to calculate 
the average and variance of most transport properties. As an example, we consider a 
quantum dot connected to source and drain reservoirs by means of two ballistic point 
contacts with A^^i and N2 propagating channels per spin direction at the Fermi level, with 
A^ = A^^i + A"2. The zero-temperature conductance is given by the Landauer formula, 
which we write as 

G = 2£!f^-t.5A5.AV (6) 




h \ N 

where S is the N x N scattering matrix and A is an A^ x A^ diagonal matrix with 

N2 if 1 < 2 < N,, 

-Ni iiNi<i<N. ^ ^ 

For large A^, the second term in equation (|^) is a small and fluctuating quantum 
correction to the classical conductance of the quantum dot. Using equations @ and 
(I), the average and variance of the conductance for A^ ^ 1 then follow as 

_2e'fN,N2 N,N2\ 

^a^^ = 7^(^J(i + W, (9) 

where the symmetry parameter /? = 1 or 2 with or without time-reversal symmetry, 
respectively. 

In the derivation of equations (P) and (P) it is important that the matrix A is 
traceless. This ensures that the non-Gaussian cumulant (^) does not contribute to 
varG, despite the fact that calculation of varG involves an average over a product of 
four scattering matrices. Similarly, the 0{N~'^) corrections to the second moment Wi of 
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equation in the presence of time-reversal symmetry do not contribute to the average 
conductance to order A^°. 

So far we have only considered elements of the scattering matrix at one value 
of the Fermi energy e (and of the magnetic field, etc.). If one wants to calculate 
averages involving scattering matrices at different energies, one needs to know the joint 
distribution of the scattering matrix iS(£:) at different values of e. To date, no full solution 
to this problem is known for > 1. However, for large N, the joint distribution of 
scattering matrix elements Sij at different values of the Fermi energy or other parameters 
continues to be well approximated by a Gaussian, while unitarity causes non-Gaussian 
corrections that are small as As before, the Gaussian part of the distribution is 

specified by its first and second moment. The first moment is zero for a quantum dot 
with ideal leads; the second moment reads || 

_ 1 ^ \ {SikSji + SilSjk) with TRS, , , 

~ AT - i(£ - £') ^ I 5ik6ji without TRS. ^ ' 

Here, and below, we measure energy in units of A/2tt, where A is the mean spacing 
between the spin-degenerate energy levels in the quantum dot without the leads. 
Equation ([T^) was originally derived using semiclassical methods ^ 



and in the Hamiltonian approach of random-matrix theory |jT2|, A derivation 

using the random scattering matrix approach is given in reference [46| and in appendix 
B. In the absence of time-reversal symmetry, the leading non-Gaussian correlations are 
described by the cumulant 

-(^n..(^i)^^*.J^'i))(^....(^2)5,;j4)) 
-(5n,.(^i)^.*.J4))(5..,.(^2)5^,J4)) 

^N-i{e,-e[)KN-i{e,-e'MN-Ke2-e[)){N-i{e,-e',)y ^ ' 
In the presence of time-reversal symmetry, 14 terms corresponding to the permutations 
^2 ^ ^2, ^1 ^ and /c2 h have to be added to equation (|TT|), respectively, as in 
equation (|^) for the energy-independent case. A derivation of equation (|TT|) is given in 
appendix B. 

Equations (0) and (|llD can be used to calculate averages and correlation functions 
for transport properties that involve scattering matrices at different energies. As 
an example, using equation (|10|) for the second moment of the scattering matrix 
distribution, the conductance autocorrelation function is found as [Hil ESl 



{G{e,)G{e,)) - {G{e,)){G{e,)) = ^^^f fef- L)^ " ^^^^ 

§ For the second moment, an exact solution was obtained using the supersymmetry approach, see 
reference [Ol. 



Scattering matrix ensemble for time- dependent transport 



7 



3. Time-dependent scattering 

For time-dependent scattering, the energies of incoming and scattered particles do not 
need to be equal. In order to describe scattering from a time-dependent scatterer, 
we use a scattering matrix S(t,t') with two time arguments. (We prefer to use the 
formulation with two time arguments instead of a formulation in which S has two 
energy arguments, since the former allows us to describe an arbitrary time-dependence 
of the perturbations.) For a quantum dot coupled to leads with, in total, scattering 
channels, the two-time scattering matrix S{t,t') relates the annihilation operators aj(t) 
and hi{t) of incoming states and outgoing states in channel i = 1, . . . , N, 

= J2 Sijit,t')sij{t')dt', 

blW= E / at(t')(5t(t',t)),,rft'. (13) 

Causality imposes that 

S{t,t') = litKt'. (14) 



Unitarity is ensured by the condition 

N 



J2 I dt{S\t" ,t)\,S,k{t,t') 

.7 = 1 



5{t" - t')5ik, 



J 

N 



I dtSij{t",t){S\t,t%k = S{t"-t')5,k, (15) 
where the Hermitean conjugate scattering matrix S^{t',t) is defined as 

{s\t',t))^^ = s;,{t,t'). (16) 

For a quantum dot without time-independent potential, the scattering matrix 
S^{t,t') depends on the difference t — t' only. (In this section, we use a superscript 
"0" to indicate that is a scattering matrix for time- independent scattering.) It is 
related the scattering matrix in energy representation by Fourier transform, 

S%t,t') = ^ deS%e)e''^'-''^/\ (17) 

Borrowing results from the previous section, we infer that the elements of S^{t,t') have 
a distribution that is almost Gaussian — the Fourier transform of a Gaussian is a 
Gaussian as well — , but with non-Gaussian correlations that are small as oo. 
Fourier transforming equation (pUj), we obtain the variance of the distribution [El] 



(5° (t, t')S',Ks, s')) = 5{t-t'-s + s')e{t - t')V\t - t') 



1 {kk5ji + ki^jk) with TRS, , , 

^ ^ 6ik6ji without TRS. ^ ^ 
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Here, time is measured in units of 27rh/ A and the function P° is given by 

V\t) = e-^^ (19) 
Fourier transform of equation ( |Tl]) gives the leading non-Gaussian contribution, 

= Kn(^iAK,,ih,t',)S'^^^^^^^^ 

X 6{ti - t[ +t2- 4 - Si + S[ -S2 + S2)9{ti - t\)e{t2 - 4)^(51 - S'l), 

with 

J'\n]T2;n;n) = [iVmin(ri,r2,r3,r4) - l]e-^(^^+^^). (21) 

Note that, in view of the delta function in equation (pOD, the function depends on 
three time variables only. Despite the redundancy, we keep the four time arguments 
for notational convenience. As before, in the presence of time-reversal symmetry, 
the expression for the cumulant is obtained by adding terms that are obtained after 
interchanging 12 ^ j2, ^ ^1, and ^2 ^ ^2, cf. Eq. (||). 

In order to calculate the defining cumulants Wi and 1^2 for the case of a chaotic 
quantum dot with a time- dependent potential, we need a statistical model for the 
scattering matrix distribution for time-dependent scattering. Such a model can be 
provided by the Hamiltonian approach P2|, or, alternatively, by extending the "stub 



model' of references to the case of time-dependent scatteringfj} In the latter 

approach, the N x N scattering matrix S(t, t') is written in terms of an M x M random 
matrix U{t, t') (with M ^ A^) and a (M — A^) x (M — A^) random Hermitean matrix iJ, 

s = pu{i- m)-^p\ = Q V^^^/^'^g. (22) 

Here P is an A^ x M matrix with Pij = 5i j and Q is an (M — A^) x M matrix with 
Qij = Si+N,j- The scattering matrices S(t,t') and U{t,t') depend on two time indices, 
and the matrix products involving U{t,t') in equation (|22| ) also imply integration over 
intermediate times. The Hermitean matrix H depends on a single time argument and 
models both time-independent and time-dependent perturbations to the Hamiltonian 
of the quantum dot. The matrix U{t,t') depends on the time difference t — t' only and 
satisfies the constraint of unitarity, equation (|T^) above. As the effect of a time-reversal 
symmetry breaking magnetic field will be included in H, cf. equation (^) below, we 
further require that the matrix U is time-reversal symmetric, 

U,,{t-t') =Uji{t-t'). (23) 

The statistical distribution of the matrix U is the same as that of the scattering matrix 
of a chaotic quantum dot coupled to a lead with M channels, but without magnetic field 

II The stub model is similar in spirit to the "quantum graph" , the spectral statistics of which is known 
to follow random-matrix theory ||50| . 
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and time-dependent potential. Hence, the first nonvanishing moments of tlie distribution 
are given by equations ( |T8| ) and ( pOD above, with replaced by U and by M. 

The physical idea behind equation ( p2D is that the time-dependent part of the 
potential is located in a "stub" (a closed lead), see figure |1|. The number of channels 
in the stub is M — A^. The matrix U is the M x M scattering matrix of the quantum 
dot without the stub; the scattering matrix S is the scattering matrix of the entire 
system consisting of the dot and the stub, taking into account the time-dependent 
scattering from the stub. The matrix R represents the time-dependent scattering matrix 
for scattering from the stub. The stub is chosen to be small compared to the quantum 
dot, so that reflection from the stub can be regarded instantaneous — that's why the 
matrix R{t) depends on a single time argument only. At the end of the calculation, 
we take the limit M — > oo. This limit ensures that the dwell time in the dot, which is 
proportional to 1/A^, is much larger than the time of ergodic exploration of the dot-stub 
system, which is proportional to 1/M. It is only in this limit that the scattering matrix 
acquires a universal distribution which is described by random matrix theory. Once 
the limit M ^ oo is taken, the spatial separation of chaotic scattering (described by 
the M X M scattering matrix U) and the interaction with the time-dependent potential 
(described by the time-dependent reflection matrix R) no longer affects the distribution 
of the scattering matrix S and the scattering matrix distribution found using the stub 
model becomes identical to that with a spatially distributed time-dependent potential 
in the Hamiltonian approach. 

A similar model has been used to describe the parametric dependence of the 
scattering matrix in the scattering matrix approach EO, H3, E^. For the parametric 



dependence of S, equivalence of the "stub" model and the Hamiltonian approach was 
shown in reference [^. The calculational advantage of the "stub" model is that, for a 
quantum dot with ideal leads, the vanishing of the first moment (Sij) = is manifest 
throughout the calculation, while it requires fine-tuning of parameters at the end of the 
calculation in the Hamiltonian approach. 

The matrix H in equation (22) can be written as a sum of three terms, describing 
three different perturbations to the Hamiltonian of the quantum dot^ 

H = V{t)l + -ffshape + -f^magn- (24) 

The first term in equation (^) represents an overall shift of the potential V{t) in the 
quantum dot. The second term represents the effect of a variation of the shape of the 

% In the Hamiltonian approach, the parameters xi{t), X2{t), V{t), and a{t) of equations (|2^)-(p7|) 
correspond to time-dependent variations of the form 

n{t) = S* + V{t)l + —^a{t)A + y {t)X,, 

where S and Xj are real symmetric random M x M matrices, j = 1, . . . , n, A is a real antisymmetric 
random M x M matrix, and 1 is the M x M unit matrix. The off-diagonal elements of these random 
matrices are Gaussian random numbers with zero mean and unit variance. The diagonal elements of 
S and Xj have twice the variance of the off-diagonal elements. 
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U(t-t') 



cavity: 

time-independent 
scattering 



R(t) 



Stub: 

' time-dependent 
scattering 



Figure 1. Cartoon of the picture behind equation (£^): Scattering from the chaotic 
quantum dot with the time-dependent potential is modeled as scattering from a 
chaotic quantum dot with time-independent potential and a stub with time-dependent 
potential. 



quantum dot, 

" X A 

i^shape(t) = Vx,(t)^ (25) 
^ — ^ TT 

i=i 

Here the Xj {j = 1, . . . ,n) are n time-dependent parameters governing the shape of the 
quantum dot, and the Xj are real symmetric random (M — A^) x (M — A^) matrices 
with ti XiXj = M'^Sij, i, j = 1, . . . ,n. Having more than one parameter to characterize 
the dot's shape is important for apphcations to quantum pumping |2T|, ^ The 



third term in equation ( p4|) represents the parametric dependence of the Hamiltonian 
on a magnetic flux $ through the quantum dot 

AA 

i^magn(t) =ia(t)^=, (26) 
TTV 2 

where A is a random antisymmetric (M — A^) x (M — N) matrix with ti A^A = M^. 
For a dot with diffusive electron motion (elastic mean free path /, dot size L ':$> I) one 
has 

where k is a constant of order unity and $ the flux through the quantum dot. One has 
K, = 47r/15 for a diffusive sphere of radius L and k = tt/2 for a diffusive disk of radius 
L |5^. For ballistic electron motion with diffusive boundary scattering, the mean free 
path I in equation ( ^Tf ) is replaced by 5L/8 and tiL/A for the cases of a sphere and 
a disk, respectively. (For the balhstic case, the value of a"^ reported in reference 
is incorrect, see reference ||5^.) In order to ensure the validity of the random matrix 
theory, the time dependence of the parameters xj and a should be slow on the scale of 
the ergodic time Terg of the quantum dot. 

Note that the description (p^-(p7D contains the dependence on a magnetic field 
explicitly. Having the full dependence on the magnetic field at our disposal, we no 
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longer need to distinguish between the cases of presence and absence of time-reversal 
symmetry. 

Expanding equation (p2D in powers of R, the scattering matrix S is calculated 
as a sum over "trajectories" that involve chaotic scattering in the quantum dot and 
reflections from the stub. Since different "trajectories" involve different channels in the 
stub at different times, each term in the expansion carries a random phase, determined 
by the random phases of the elements of U. Hence, elements Sij will have a distribution 
that is almost Gaussian for large A^, since they are sums over many contributions with 
random phases. Unitarity, imposed by the constraint (|15|) for the matrix U and the 



form of the matrices S and R in equation (|2^) R, leads to corrections to the Gaussian 



distribution that are small as becomes large. 

The Gaussian part of the distribution of the time-dependent scattering matrix 
S{t,t') is specified by the second moment, 

Wl^-''\tAs,s') = {S,^{t,t')S;,{s,s')), (28) 

whereas the leading non-Gaussian corrections are described by the cumulant 

VV2 [ti,ti,t2,l2: Si^ S-^^, S2, S2) 

= ('5nii(tl,t'l)5i2i2(t2,4)'5fc,,,(Sl,s'i)54;^(s2,4)) 

- {<Shh{tl,t[)Sl^i^{Si, s[)){Si^j^{t2,t'2)Sl^l^{s2, S2)) 

-(5.,,,(ti,t;)5,*,j.2,4))('5.2.2(i2,4)5,*^j5i,/i)). (29) 

The central result of this paper is a calculation of the cumulants Wi and W2 for 
time-dependent scattering. Details of the calculation are reported in appendix C. For 
the second moment Wi we find 

Wl'-'^'it, t'; s, s') = 6{t -t'-s + s')e{t - t') [6,u6,{D{t, t'- s, s') + 6a6jkV{t, t'; s\ s)] , (30) 
with 

Pfr, a; r', a') = exp 



N\T-a\-'^-^ jy di {y{a + r,0 - V{a' + ^0) (31) 



and ?7 = sign {r — a), rj' = sign (r' — a'). The first term in equation (PI1| ) is the analogue 
of the diffuson from standard diagrammatic perturbation theory, while the second term 
corresponds to the cooperon. For notational convenience, both terms are denoted by 
the same symbol V. (Note that the order of the time arguments s and s' is reversed in 
the second term of equation (|30|) .) The leading non- Gaussian corrections are given by 
the cumulant W2 for which we find 

w^'''^''''-'''''''''\t,,t[-, t2, t'2; Si, s[; S2, s'2) = e{t, - t\)e{t2 - t'2)e{s, - s[)e{s2 - 4) 

X ^hk^^hl2^i2k2^j2lMh - t'l +t2- t'2 - Si + s[ - S2 + S2) 

X J^(ti,t'i; t2, t'2; Si, s'l; S2, s'2) + permutations. (32) 
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S2 



,0 — <i 



X2 




S2 t, , S2 t] 
t'2 Si S'l 





'^t'2S^ {^2^2^, 
t2 si ^t'2 Sj 

S2 t, ^S2 t, 

sit'ij \sit\ 

t2S?\ Ft2Sl^, 




Figure 2. Diagrammatic representation of the contributions to the leading non- 
Gaussian correlator W2 of equation (p9[), which involves four scattering matrices. 
The top left diagram has weight J^{ti,t'i;t2,t'2; si, s'l, S2, s'2), cf. equation (^2|). The 
top right diagram is obtained by interchange of two scattering matrices and has 
weight T(ti,t'i;t2,t'2; S2, s'2; si, s'l). These two diagrams give all contributions to the 
cumulant W2 in the absence of time-reversal symmetry. In the presence of time-reversal 
symmetry, the fourteen lower diagrams, corresponding to the reversal of one or more 
directions of the vertices, contribute as well. 



The "permutations" in equation (|3^) refer to one term corresponding to the permutation 
{si, s[, ki,li) ^ (s2, S2, A;2, ^2) of the third and fourth arguments of W2 and fourteen 
more terms corresponding to the interchange of incoming and outgoing channel and 
time arguments within the second, third, and fourth argument of W2- A diagrammatic 
representation of the cumulant (^) and the relevant perturbations is shown in figure 
The kernel J-' reads 

J^in, cTi; T2, CT2; t[, a[; t^, a'^) = j d^V{ai + r]i^, a^, a[ + r][^, a[) 

X V{t2 -T[ + cr[ + r]2^, 02] - n -f ai -f 7^2^, ^2) 
X I?(ri, (Ti ?7i^; T'^-Ti + ai+ r]^^) 
X ViT2, T2-T[ + a[+ ri2^; t[, a[ + ri[0 



X 



N + 4 ^[s„(cTi + r/iO - Xm{r2 - n + ai + r]'^^)] 

m 



X [Xmicr'i + r][^) - x„(r2 - t[ + cr[ + r]20] - ^(^ - In " ^^il) 
- - \t[ -a[\) + 2[7],a{a, + r/iO " ^2«(^2 - n - + 
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X [r][aia[ + r][^) - r^^air^ -T[ + a[+ r^^O] | , (33) 

where we abbreviated rji = sign (ri — cri), 772 = sign (r2 — 0-2), ri[ = sign (r{ — a[), 
and ri'2 = sign (rg — (T'2). Note that equations (p0|)-(|33|) cover both the cases with and 
without time-reversal symmetry through the exphcit dependence on the magnetic flux a. 
If time-reversal symmetry is fully broken, all permutations in equation (^2]) that involve 
the interchange of incoming and outgoing channels, corresponding to the fourteen lower 
diagrams in figure 0, vanish, and only the first two diagrams in figure ^ remain. Partial 
integration of the intermediate time ^ allows one to rewrite terms between brackets {. . .} 
in equation (|33D, see appendix C for details. Finally, one verifies that the result (^) is 
recovered for a = and a ^ 1 , corresponding to presence and absence of time- reversal 
symmetry, when the parameters Xj and a do not depend on time. 



4. Applications 

In order to illustrate the use of equations (|30D-(p3D, we return to the example of section 
and consider transport through a chaotic quantum dot coupled to two electrons 
reservoirs by means of ballistic point contacts with A^^i and N2 channels, respectively. 
The scattering matrix of the quantum dot has dimension = A'^i + A^2- The current 
through the dot is defined as a linear combination of the of the currents through the 
two point contacts, 

N 

m = e J2 {4m^Mt) - bl(t)A.,b,(t)) , (34) 

where the N x N matrix A was defined in equation (|^) and the operators aj(t) and bj(t) 
are annihilation operators for incoming and outgoing states in channel i = 1, . . . , N in 



the leads, respectively, see section ^. The advantage of the definition (|3^) for the current 
through the quantum dot, instead of a definition where the current through one of the 
contacts is used, is that it simplifies the ensemble average taken below. Both definitions 
of the current give the same result for the quantity of interest, the integral of I{t) over 
a large time interval t\ < t < tf. 

The electron distribution function for the electrons entering the quantum dot 
from the leads is given by the Fourier transform f{t) of the Fermi function in the 



corresponding electron reservoir ||57|, E8l 59 



a](f)a,(t) =/.,(t'-t), 

a,(f)al(t) = 4-(t-t'), (35) 
de e^^l^ 



where we defined 



fij (^) ^ij 

= s. 



2hsmh{TTkTt/h) 
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Ut) = 6,,6{t)-Ut). (36) 

Here /Xj is the chemical potential of reservoir 1 for 1 < i < A^^i and the chemical potential 
of reservoir 2 for Ni < i < N. 

Substitution of equations (^) and (|13|) into equation ( PH ) allows us to calculate 
the time-averaged expectation value of the current through the quantum dot for a time 
interval ti < t < tf, 

I = r dt ! dtidt2 tr [5{t - ti)M{t - ta) - S\ti, t)AS{t, ts)] f{ti - ta). (37) 

- ti Jt. J 

(A factor two has been added to account for spin degeneracy. The time interval 
ti < t < t{ during which charge is measured is taken to be the largest time scale in 
the problem.) 

In the absence of a source-drain voltage, equation (0) describes the current that 
is "pumped" by the time-dependent potential in the dot, 

2e i'^' 



'pump 



ti - t 



- dt J dtidt2tTS\ti,t)AS{t,t2)feci{tl-t2), (38) 



where /eq is the Fourier transform of the Fermi function. Equation (38) was first derived 



in reference 35 ; it reduces to the current formulae of references [531 and [541 in the 



adiabatic limit, where the time dependence of the potential of the quantum dot is slow 
compared to the dwell time in the quantum dot. At small bias voltage, there is a current 
proportional to the bias, / = GV, where G is the (time-averaged) conductance of the 
dot. The conductance G can be calculated from equation ( pTj) by setting fii = AacV 
and then linearizing in V 



G 



2e^ 
h 



mN2 



27ri 



- dt j dtidt2{ti-t2)tTAS{t,ti)AS^{t2,t)f^^{ti-t2 



N -^^y -^^.-.v^. ^..j -(39) 

Here /eq is the Fermi function in the absence of the external bias. For time-independent 
transport, equation (|39|) is equal to the Landauer formula (|^). 

Conductance. The ensemble average and the variance of the conductance G for 
a quantum dot with a shape depending on a single time-dependent parameter x was 
calculated by Vavilov and Aleiner using the Hamiltonian approach [^. Using the 
scattering matrix correlator (0), their result for (G) is easily reproduced and generalized 
to arbitrary values of the (time-independent) magnetic field. 



2e'N,N2 
(G) = -jj^ + SG, 

2e^NiN2 
~ hN{ti - ti) 

X exp 



(40) 



6G 



dt 



dr 



{N + W)t -2 dTi{x{t - r + Ti) - x{t - ti))' 



(41) 



The correction term 5G of equation (|41| ) is the weak localization correction; it results 
from the constructive interference of time-reversed trajectories. The presence of a time- 
dependent potential breaks time-reversal symmetry and suppresses the weak localization 
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correction. Vavilov and Aleiner investigated the case x{t) = Sxcos{ujt) of a harmonic 
time dependence for the parameter x in detail In that case, the suppression of weak 
locahzation increases with increasing frequencies and saturates at a value 

if {6xf < + Aa\ 



6G = — X 



N'h ] ./^ if(5x)2»A^ + 4a^ 



(42) 



for frequencies hw > A^A (Apphcabihty of random matrix theory requires that 

to <^ 1/Terg, where Terg is the time for ergodic exploration of the quantum dot.) If the 
fluctuations of the parameter x are fast and random on the scale 2'irh/NA of the delay 
time in the dot, they may be considered Gaussian white noise, 

{x{t)x{t')) = \^6{t-t'). (43) 



In that case, the exponent in equation (^) can be averaged separately, and one finds 
the result 

SG = ?^!iM^. (44) 

The same suppression of weak localization was obtained previously to describe the 
decohering effect of the coupling to an external bath |HT|, Note that the 

strong-perturbation asymptote for white noise is different from the strong-perturbation 
asymptote for fast harmonic variations of the dot's shape. The cause for this difference 
is the existence of small time windows in which time-reversal symmetry is not violated 
near times t with cos{ujt) = ±1 for harmonic variations oc cos{ujt), while for a random 
time dependence of x{t) no such special times around which time-reversal symmetry is 



preserved exist p5|, P6| . 

Similarly the variance of the conductance can also be expressed in terms of the 
correlator (|28|). (As in the time- independent case, the non-Gaussian correlator (^) 
does not contribute to the variance of the conductance.) We refer to reference 
for the detailed expression for varG and an analysis of the effect of a harmonic time 
dependence of the shape function x{t). Conductance fluctuations for the case when x{t) 
is a sum of two harmonics with different frequencies were considered by Kravtsov and 
Wang 

Pumped current. To first order in the pumping frequency u, the current in the 
absence of a source-drain voltage is nonzero only if two or more parameters Xj that 
determine the dot's shape are varied independently. Even then, the ensemble average of 
the pumped current is zero, and the first nonzero moment is (/^). The ensemble average 
(P) was calculated in reference for small pumping amphtudes xi(t) = Sxi sin(a;t), 
X2 = 5x2 sin(tt;t + 0), 

/r2\i/2 euj5xi5x2 . , 

^I)' = 2nN ^^"'^' ^^^^ 
independent of the presence or absence of a magnetic field. The case of pumping 
amplitudes of arbitrary strength was considered in reference ^4 . 
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Figure 3. Diagrams representing the two contributions to the variance of the pumped 
current. Following the notation of reference js^, dotted lines correspond to the 
scattering matrix S, thick solid lines correspond to the fixed matrices A and /, and 
thin solid lines correspond to the Kronecker delta's in the average over the ensemble 
of scattering matrices. 



Beyond the adiabatic regime, one time-dependent parameter is sufficient to generate 
a finite current through the dot The second moment {P) in that most general 



case was first calculated in reference , using the random Hamiltonian approach. The 
second moment, which involves an average over four scattering matrix elements, can 
also be obtained in the random scattering matrix approach, using equations (pOD-(p^) 



of the previous section. We then find that there are two contributions to (/pump)- 
contribution with two Gaussian contractions of scattering matrices (giving a factor Wf) 
and one contribution which involves a correlator of four scattering matrices (giving a 
factor W2)- Diagrams representing these two contributions are shown in figure |^. Adding 
both contributions, we find 



2 8e^iViiV2 

\ pump; ^^^^ _ ^.p 



ti JO 

X V{t,t - T - T';t',t' - T 



dtdt' / drd^d^' / rfr7eq(2r')/cq(-2r') 



X V{t - 
X V{t' 



T 



T')V{t', t' -T + T';t,t-T + t') 
i;t-T + T',t-T + T' -i) 



T + r'A' 



t'A' 



T 



T 



{6(0 - N){6{e) - iV) + 4 J2Mt' -r + r')- x^t' 



X {xmit -r + r) - Xmit - r - r')) 



-')) 



(46) 



independent of the value of the magnetic field. Using equation ( |(J.4D of appendix C to 
express the delta functions in terms of the functions Xm{t) and a total derivative of V, 
performing partial integrations, and shifting t^t — r,t'—>-t' — T, this can be rewritten 
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as 



(J2 ) 

\ pump/ 



32e^NiN2 
N{t( - ti)2 7o 



tf-T 



drdidi' / rfr' / rftrft7eq(2r )/eq(-2r') 



X ©(t + r, t - r'; t' + T,t' - T')V{t - t' ,t - t' - ^^t + t' ,t + t' - ^ 
X V{t' + T,t' + T';t + T,t + T')V{t' + r', t' + r' - t' - r', t' - r' - 



X 



m,n \ ± J \ ± ) 



+ 5^ 5^ ±a;™(t' - r ± r') ±a;™(t - r ± r') 



(47) 



This expression agrees with the resuh found by Vavilov, Ambegaokar, and Aleiner pS 
We refer to reference for a detailed analysis of equation (^) for the limiting cases 
of adiabatic pumping and high-frequency pumping with one and two time-dependent 
parameters. 

Noise. The current noise is defined as the variance of the charge transmitted 
through the quantum dot in the time interval U <t <ti 

1 



S 



tf-ti 



dtdt' {I{t)l{t') -I{t)I{t'] 



(4^ 



As in equation (0), denotes a quantummechanical or thermal average, not an 
ensemble average. Performing the quantummechanical and thermal average over the 
incoming states [0, |5^, |59|, the noise power S can be calculated as 
2e2 ■ 



S 



tf-ti 



dtdt' j dtidt2dt\dt'^ii t)A5(t, ta) - - t)M{t - ts)) 

X f{t\ - t^) (5t(t;, t')KS{t\ t'^) - - t')M{t' - t',)) f{t, - t',] 



(A factor two has been added to account for spin degeneracy.) In the absence of a time- 
dependent potential. Equation P9| ) represents the sum of Nyquist noise and shot noise 
. With time-dependence, it contains an extra contribution to the noise that is caused 
by the time dependence of the potential in the quantum dot ^ ^ 

Averaging equation (^) for an ensemble of chaotic quantum dots, we find 

{S) = 3'' + S^ + 5^, 

= 2kTh{G), 



= eVh{G) 



coth 



eV 2kT\ 



27im V 2kT eV J ' 



S 



p _ e^NiN2{kT/hf 
2N{t( - to 



dtdt' 



1 



V{t,t-i-t\t' -i)di 



X 



- 2{Nl + A^l) ?,in^[eV{t - t') /2h] 
sinh^[7rA;T(t -f)/^] ' 



(49) 



where {G) is the average (time-dependent) conductance, see equation (^0|) , and V = 
(/ii — /i2)/e the bias voltage. The above ensemble averages for the Nyquist noise and shot 
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noise are the same as the noise power found in the absence of a time-dependent potential 
||70| , up to an eventual weak localization correction. The extra noise generated by the 
time-dependence of the dot shape is fully described by the term S^. In the adiabatic 
regime hw <^ A^A, the pumping noise can be written as 



e^NiN2{kT/hf 
2N{t( - ti) 



dt'dt 



1 

iV2 



X 



- 2{Nf + Nl) sm^ [eV{t-t')/2h] 
'fj/h] ■ 



9r . . , . , (50) 

In the absence of a bias voltage, eV = 0, Equation (|50|) has been analyzed in detail 
by Vavilov and the authors in reference [0. For one time-dependent parameter 
x{t) = 6xcos{ujt) it is found that 



S 



(51) 



ue^N^N2 \ 27i{5xy (coth ^ - ^) , if (Sx)^ < max(l, k^T^/h^u^) 
Ti^N^ [ 3\Sx\N'/\ if {6xy > Armax(l, eT^/h^u^). 

An applied bias voltage has an effect on the pumping noise only if eV > 
max(/ia', fcT, /ia;|5x|/A^^/^). And even then, the effect of the applied bias is limited 
to a reduction of 5''^ by a numerical factor 2NiN2/N'^. In this respect, the effect of an 
external bias on the pumping noise is much weaker than that of temperature, which 
tends to suppress as soon as kT > hiu max(l, \5x\/N^^'^) 



37 



5. Conclusion 



In summary, in this paper we have extended the scattering approach of the random- 
matrix theory of quantum transport to the case of scattering from a chaotic quantum dot 
with a time-dependent potential. We addressed the limit that the number of channels 
A^ coupling the dot to the electron reservoirs is large. In this limit, the elements of 
the scattering matrix have a distribution that is almost Gaussian, with non-Gaussian 
corrections that are small as A^ becomes large. We calculated the second moment, 
which defines the Gaussian part of the distribution, and the fourth cumulant, which 
characterizes the leading non-Gaussian corrections. 

The advantage of the scattering matrix approach is that, once the scattering matrix 
distribution is calculated, the computation of transport properties is a matter of mere 
quadrature. As an example, we calculated the conductance of a quantum dot with a 
time-dependent potential or the current pumped through the dot in the absence of an 
external bias, and found agreement with previous calculations of Vavilov et al. that 
were based on the Hamiltonian approach |2^, The results derived here were 

used for the calculation of the current noise generated by the time- dependence of the 
potential in the quantum dot by Vavilov and the authors 1^^. The current noise in the 
presence of both a time-dependent potential in the dot and a bias voltage was studied 
here. 

Whereas the first four moments of the scattering matrix distribution that we 
calculated here are sufficient for the calculation of most transport properties — most 
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transport properties are quadratic or quartic in the scattering matrix — , we need to 
point out that there are observables that cannot be calculated with the results presented 
here. First, in the presence of one or more superconducting contacts, (averaged) 
transport properties may still depend on higher cumulants of the distribution, despite 



the fact that these are small by additional factors of [jTT|. Second, the results 



presented here fail to quantitatively describe transport properties for very small A^, 
which can have strongly non-Gaussian distributions. Further research in these directions 
is necessary. 
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Appendix A. Nonideal contacts 

Nonideal contacts are characterized by channels that have a transmission coefficient 
Tj smaller than unity, j = 1,...,N. The imperfect transmission of the contacts is 
characterized by an x reflection matrix r^t, t'), for which we take the simple form 

r,{t,t') = {l-rY/'6{t-t'), (A.l) 

where F is an x diagonal matrix containing the transmission coefficients Tj on 
the diagonal. The direct backscattering from the contacts is fast compared to the 
scattering that involves ergodic exploration of the dot, hence the delta function 6{t — t') 
in equation ( |A.1| ). In order to describe time dependent scattering with nonideal leads. 



we use a modification of the stub model of equation (^) |]7T], pn| , 

5 = r, + ri/2^flri/2^ s"" = pu{i- m)-^p\ (a.2) 

The first term in equation ([A.2|) takes into account the direct backscattering at the 
contact for electrons coming in from the reservoirs, whereas the extra term in equation 
( |A.3| ) describes backscattering at the contact for electrons coming from the dot. The 
additional factors F^/^ in the second term of equation ( |A.2| ) account for the decreased 



transmission probability for entering or exiting the quantum dot. With the inclusion 
of refiection in the contacts as in equation ( [A.2|) , the scattering matrix approach for 
time-independent scattering was proven to be fully equivalent to the Hamiltonian 



approach with arbitrary coupling to the leads [|20[. The corresponding distribution of 
the scattering matrix S for time-independent scattering is known as the Poisson kernel 

Like in the case of ideal leads, the distribution of the elements of the scattering 
matrix S for a quantum dot with nonideal leads is almost Gaussian, with non-Gaussian 
corrections that are small if A^ ^ 1. The main difference with the case of an ideal 
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contact is that, as a result of the direct reflection from the contact, the average of S is 
nonzero for a nonideal contact. The fluctuations of S around the average are described 
by <S^, cf. equation ( |A.2| ). In order to find the distribution of S^, we note that the 
expression ( |A.2| ) for is formally equivalent to the stub model equation ( p2D used to 
describe time-dependent scattering from a quantum dot with ideal contacts. Hence we 
conclude that the moments of can be obtained directly from the results for the case of 
ideal contacts, see section | and appendix C, provided we substitute equation ( |A.3[ ) for 
the matrix R. This amounts to the replacement 5 ^ 5^ in the final results (^)-(|33|), 
Tj in equation (^l]) , and r| in equation (|33D . 

Appendix B. Correlators for time-independent scattering 

The scattering matrix correlators for time-independent scattering serve as input for the 
calculation of the correlators for time-dependent scattering. They can be calculated 
using the Hamiltonian approach (see references fl^ , |4^ , |45[|), or, alternatively, in the 
scattering matrix approach, using a time-independent version of the "stub model" of 
section ^. Following the latter method, the scattering matrix S is written as 

S{e) = PU{1 - KU)-^P\ R = Q^e^'^'^/^^^Q. (B.l) 

Here the matrices P and Q are as in equation (|22|), whereas W is an M x M unitary 
matrix taken from the circular orthogonal ensemble or circular unitary ensemble of 
random matrix theory, depending on the presence or absence of time-reversal symmetry. 
The picture underlying equation ( [B.l| ) is that a stub with M — N scattering channels is 
attached to the chaotic quantum dot as in figure ^ such that the dwell time in the stub 
is much larger than the dwell time in the dot, but much smaller than the total dwell 
time in the combined dot-stub system. The first condition implies that the M x M 
scattering matrix of the chaotic dot (without stub) may be taken energy independent, 
and distributed according to the appropriate circular ensemble from random matrix 
theory. The total scattering matrix S then acquires its energy dependence through the 
energy dependence of the (M—N) x (M—N) reflection matrix R of the stub. The second 
condition, which requires M ^ N, ensures that the dot plus stub system is explored 
ergodically before an electron escapes into the lead, so that the spatial separation of the 
energy dependence (stub) and chaotic scattering (dot)) does not affect the correlators 
of the scattering matrix 5.0 



Using the diagrammatic technique of reference [36| to average over the random 



+ Note that this version of the "stub model" is different from that used in the main text. In time 
representation, the matrix 14 of equation ( [B.l| ) is proportional to a delta function S{t — i'), whereas 
the matrix R involves a time delay with time t — t' = 2Trh/MA. For the model of section || of the 
main text, the time delay is described by U, whereas scattering from the stub is instantaneous. Both 
versions of the "stub model" are equivalent to the Hamiltonian approach. Which one to use is a matter 
of convenience. 
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unitary matrix U, we find that the second moment Wi is given by 

urirM(^. /X ^ I X / i^ik^ji + ^ii^jk) with TRS, 

^ ^ ' ^ M -tTR{e)R^e') \ 6ikSji without TRS. ^ 

Substitution of equation ( |B.1| ) for R gives equation ( |10|) of section |[ Note that equation 
(p!0|) is vahd in the semiclassical hmit of large only. Within the diagrammatic 
technique this follows from the observation that for large the only contributions 
to Wi are the "ladder" and "maximally crossed" diagrams, whereas for small 
more contributions exist and a non-perturbative calculation is needed to calculate 
the scattering matrix correlator |]5B|. The correlator Wi{e,e') was calculated by 
Verbaarschot et al. in reference |T2| for arbitrary using the Hamiltonian approach 
and the supersymmetry technique. 

For the cumulant W2 we find in the absence of time-reversal symmetry 

X [M -tiR{ei)R{e2)R\e\)R\e'2)] 

X [M - ti R{ei)R{e\)]-^[M - ii R{e2)R{e\)]-^ 

X [M - tri?(ei)i?(e^)]-^[M - tr i?(e2)i?(4)]"^ (B.3) 

In the presence of time-reversal symmetry, fourteen terms corresponding to the 
permutations 12 j2, ki ^ li, and k2 ^ I2 have to be added. Equation (|lT]) is 



recovered upon substitution of equation ( |B.1|) for R 



Appendix C. Correlators for time-dependent scattering 



In this appendix we present the derivations of equations (|30D-(p3D. 

We first calculate the second moment Wi of the scattering matrix distribution, 
equations ( pOD and (|3lD. To find Wi we use equation (|2^ ) to expand S in powers of 
U and R and then average over U. In the limit of large M and large A^, that average 
can be done using the cumulants ([T8|) and (|20|) and the diagrammatic rules of reference 
36| . This calculation is similar to the standard diagrammatic perturbation theory: the 



matrices U, U\ and R{t) play the role of the unperturbed retarded and advanced Green 
functions and the random potential, respectively. 

Performing the average over lA this way, we find that, to leading order in and 
A^~^, the cumulant Wi is dominated by two leading contributions: the "ladder diagram" 
of figure Cl and the "maximally crossed diagram" of figure 02. Since every factor in 



these two diagrams involves equal time differences for S and 5*, we conclude that this 
contribution to Wi is nonzero only ii t — t' = s — s', cf. equation ([T8|). Further, we 
conclude that the ladder diagram gives a nonzero contribution only if i = k and j = I, 
while the maximally crossed diagram contributes when i = I and j = k. 

We first consider the contribution of the ladder diagram, which we write as 

6ik6ji6it -t'-s + s')D{t, t'- s, (C.l) 
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Figure CI. (a) Notations, following reference (b) Calculation of the kernel 

D{t, t'; s, s'). (c) Diagrammatic representation of the Dyson equation ( |C.2|) for D{t,t — 
t; s, s — t), with t' = t — T,ti=t — ti, s' = s — t, and si = s — ti. (d) Diagrammatic 



representation of the differential equation (C^) for D{t, t — t; s, s — t). Left hand side: 
{M +dr)D{t, t-r; s, s-r); Right hand side: 5(T)+tr R{t~T)W {s-T)D{t, t-T\ s, s-t). 



where the kernel D is the equivalent of the "diffuson" from standard diagrammatic 
perturbation theory. Note that, in view of the delta function in equation ( |C.1| ), the 
kernel D depends on three arguments, not on four. For notational convenience, we 
prefer, however, to continue to use the two initial times t' and s' and the two final times 
t and s to denote the time- arguments of D. 

Considering the ladder diagrams to all orders, the diffuson D is found to obey the 
Dyson equation 



(C.2) 



D{t, t~T;s,s-T) = e{T)e'"'^ + ^(r) / dnD{t, t - n; s, s - ri 

Jo 

X ti R{t - Ti)R^{s - ri)e-^^(^-^i). 
The solution of equation (|C.2D is 



D{t' + r, t'; s' + r, s') = ^(r) exp 



dn {M - tr R{t' + n)R\s' + ri)) 



(C.3) 



where we used that D = if t < 0. Substitution of i? = exp(27riJ/A) reproduces the 
first term in the result (|30|). For future use, we note that the function D of equation 
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+... 



Figure C2. Diagram representing the kernel C(t, t'] s, s') of equation (C.7) 



( |U.3| ) obeys the differential equations 

^D{t, t-T;s,s-T) = 5{t) -[M-ti R{t - t)R\s - r)] D{t, t-T;s,s-T), (C.4) 



-^D(t' + r, t'; s' + r, s') = 5(t) -\M-tr Rit' + t)RHs' + r)l Dit' + r, t'] s' + r, s'). 

(C.5) 

Calculation of the contribution of the maximally crossed diagram proceeds in an 
analogous way. This contribution reads 

6{t-t'-s + s')6u6,kC{t, t'- s, s'), (C.6) 

where the analogue of the Cooperon is given by 



C{t' + T,t';s' + T, s') =^(r)exp 



dn (M - ti R{t' + Ti)R\s' + T- ri)) 



(C.7) 



Substitution of i? = exp(27r/f/A) gives the second term of equation 

We now turn to the four scattering-matrix correlator (^), which is the equivalent 
of the Hikami box in standard diagrammatic perturbation theory. We first calculate the 
first term of equation (^2]). It is represented diagrammatically in figure |C3| . There 
are two contributions: One contribution involving Gaussian contractions with the 
cumulant (^) only, which is depicted as the first term on the r.h.s. of figure |C3| , and 
one contribution that involves the non-Gaussian contraction of equation ( ^Of ) once and 
otherwise Gaussian contractions, see the second term on the r.h.s. of figure |C3| . Adding 
those two contributions, the function is found to be, 

J^{ti,t[; t2, t'^; Si, s[; S2, s'^) = J dTD{t[ + r, t[; s[ + r, s[) 

X D{t2 - si + s[ + r, 4; S2-ti+ t[ + r, 4) 

X D{ti, t[ + r; S2, S2 - ti + t[ + T)D{t2, t2 - Si + s[ + r; Si, s[ + r) 
X ti R{t[ + t)R\s2 - tl + t[ + T)R{t2 -si + s\ + t)R\s\ + r) 

- j dTidT2dT3dTi{M + dr^)D{t\ + ri, 4; s\ + ri, s\) 

X (M + dr,)D{t'2 + r3,4; ^ + rg, 4)(M + 9,jZ^(ti,ti - Ts; S2, S2 - T2) 
X (M + dr^)D{t2, t2 - si, si - T4) 

X 9{ti - t[ - Tl- T2)9(t2 - 4 - ^3 - r4^)0{si - s[ - Tl- T4)9{s2 - 4 - T2 - Tg) 

X J^^iti - t[ - Tl- T2] t2-t'2-T3- T4; Si - s'l - Ti - T4] S2 - 4 - T2 - T3) . (C.8) 
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Si 



«2 



: = 



t2 




Si 



Si 



Figure C3. Diagrams representing the two contributions to the correlator 
J^(ti,t'i;t2,t'2; si, s[; S2, s'2) (l.h.s.). The first diagram on the r.h.s. contains Gaussian 
contractions only; the shaded blocks in denote the kernel D. The second diagram on 
the r.h.s. contains one non-Gaussian contraction, which is in the center of the diagram; 
the shaded blocks represent factors of the form (M + dr)D, see equation ( |C.8[) and 
figure Cli. 



Here we used equation (|C.5|) to express the four legs of the diagrams of figure |C3| b in 
terms of the diffuson D and its derivative. The second term in equation (|C.8|) can be 
simphfied noting that the time integration is dominated by all four arguments of 
being of order 1/M. Using the smallness of these time arguments, the diffusons can be 
expanded around ti—t']^ — ri — r2 = — ^2^T3 — T4 = Si — s^ — ti — t^ = S2 — S2 — T2 — T3 = 
and three of four time integrations can be done. The result is 

J^(ti, t[; t2, t'21 Si, s'^] S2, S2) = J dTT(ti,t[; t2, t'2] Si, s\] S2, T)D{t[ + r, t[; s[ + r, s[) 

X D{ti, t[ + r; S2, S2 - ti + t[ + T)D{t2, t2 - Si + s[ + r; , Si, s[ + r) 

X D{t2 -si + s[ + T, t'2; S2-ti+ t[ + r, 4), (C.9) 

where we abbreviated 

r = M - ti R{t[ + t)R^{s2 - ti + t[ + r) - tr i?(t2 - Si + s[ + t)R\s[ + r) 
+ ti R{t[ + r)R^{s2 - ti + t[ + T)R{t2 - si + s[ + t)R\s[ + r) 
-5{t^-t\-T)-5{si-s\-T). (C.IO) 

We used equations ( |C5| ) and ( p.4| ) to calculate time derivatives of the diffusons. 
Alternatively, using a partial integration, the function can be expressed by equation 
(^) with 

r = M - ii R{t\ + t)R\s\ + r) - ii R{t2 -si + s\ + t)R\s2 - ti + t[ + r) 
+ ti R{t\ + t)R\s2 -ti+ t[ + T)R{t2 -si + s[ + t)R\s[ + r) 



- 5{t) - 5{t -ti+ t[ + S2 



(c.ii; 



or with r given by a convenient linear combination of equations (IClOj ) and ( C.ll|) with 
coefficients Ci, C2 satisfying the condition Ci + C2 = 1. 

Finally, using equation (|22|) for R, the first term of equation (|32D is obtained. The 
other contributions to equation (32) can be found after permutation of the channel 
indices and time- variables as indicated in figure ^ 
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